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We report on a search for charged massive resonances decaying to top (t) and bottom {b) quarks in 
the full data set of proton-antiproton collisions at center-of-mass energy of ^/s = 1.96 TeV collected 
by the CDF II detector at the Tevatron, corresponding to an integrated luminosity of 9.5 fb“^. No 
significant excess above the standard model (SM) background prediction is observed. We set 95% 
Bayesian credibility mass-dependent upper limits on the heavy charged particle production cross 
section times branching ratio to tb. Using a SM extension with a W' —>■ tb and left-right-symmetric 
couplings as a benchmark model, we constrain the W' —>■ tb mass and couplings in the 300 to 900 
GeV/c^ range. The limits presented here are the most stringent for a charged resonance with mass 
in the range 300 - 600 GeV/c^ decaying to top and bottom quarks. 


Several modifications of the standard model (SM) of 
particle physics predict the existence of massive, short¬ 
lived states decaying to pairs of SM leptons or quarks. 
Such a resonance decaying to a top (t) and a bottom (5) 
quark, tb, appears in models such as left-right-symmetric 
SM extensions [T], Kaluza-Klein extra dimensions ma, 
technicolor n la or little Higgs scenarios featuring 
one or more massive charged vector bosons, generically 
denoted as W. Searches for W' bosons in the W -T 
tb decay channel are complementary to searches in the 
leptonic decay channel W -T Iv, and probe the most 
general scenario where the couplings of the W boson to 
fermions are free parameters. 

Recent searches in the W -T tb channel have been 
performed by the CDF [7] and DO [5] collaborations at 
proton-antiproton {pp) collisions at 1.96 TeV center-of- 
mass energy (CM) at the Tevatron, and by the ATLAS 
[9] and CMS [10] collaborations in proton-proton collisions 
at 8TeV CM energy at the Large Hadron Collider (LHC). 
For mass scales approaching and surpassing 1 TeV, the 
LHC experiments have superior sensitivity to the Tevatron 
experiments due to the enhancement of the production 
cross section at the higher center-of-mass energy of the 
collisions. However, in the mass region well below 1 TeV 
the Tevatron experiments have greater sensitivity due to 
the relative suppression of gluon-initiated backgrounds 
compared to the quark-initiated signals such as the one 
under consideration here. 

In this Letter, we present a novel search for charged 
massive resonances decaying to tb quark pair. The search 
is performed in events where the top quark decays to a 
Wb pair and the W boson decays to a charged lepton and 
a neutrino; the two bottom quarks hadronize and produce 
two clusters of particles (jets). Since no assumptions on 
the signal model other than on the natural width are 
made, this search is sensitive to any narrow resonant 


state decaying to a tb final state. A simple left-right 
symmetric SM extension m, predicting the existence of 
W bosons of unknown mass and universal weak-coupling 
strength to SM fermions, is used as a benchmark model. 
The reconstructed width of the signal is dominated by 
resolution effects; the test signal is therefore applicable 
for any FF'-like particle whose width is small compared 
to the experimental resolution. 

The data were collected at the Tevatron pp collider at 
a center-of-mass energy of 1.96 TeV and were recorded 
by the CDF H detector (T^j. The detector consists of a 
silicon microstrip vertex detector and a cylindrical drift 
chamber immersed in a 1.4 T magnetic field for vertex 
and charged-particle trajectory (track) reconstruction, 
surrounded by pointing-tower-geometry electromagnetic 
and hadronic calorimeters for energy measurement, and 
muon detectors outside the calorimeters m- 

We analyze events accepted by the online event se¬ 
lection (trigger) that requires either the event missing 
transverse-energy to satisfy 45 GeV or, alterna¬ 
tively $x > 35 GeV and the presence of two or more jets, 
each with transverse energy Ep > 15 GeV. The full data 
set corresponds to an integrated luminosity of 9.5 fb“^. 
Offline, we select events with > 50 GeV, after cor¬ 
recting measured jet energies for instrumental effects [T4|. 
We further require events to have two or three high-A-p 
jets, where the two jets ji,j 2 with the largest transverse 
energies, and E^ , are required to satisfy E^ > 35 
GeV and E^ > 25 GeV; the jet energies are determined 
from calorimeter deposits and corrected using charged- 
particle momentum measurements m One leading jet 
is required to be within the silicon detector acceptance, 
|? 7 | < 0.8; the other satisfies |? 7 | < 2.0. In addition to 
the large missing transverse-energy indicating the pres¬ 
ence of a high-pT neutrino, the presence of a IF boson 
decaying to an ez/g or pair is confirmed by requiring 
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a reconstructed electron or muon. Leptonically decaying 
r leptons are collected in the same way. Hadronically 
decaying r leptons from the W decay chain are mostly 
reconstructed as jets in the calorimeter. Three-jet events 
are thus retained, while events with more than three jets 
with > 15 GeV and |? 7 | < 2.4 are excluded. The 
majority of the background at this stage is quantum crho- 
modynamics (QCD) production of multijet events, which 
yields generated through jet-energy mismeasurements. 
Neutrinos produced in semileptonic 6-hadron decays also 
contribute to the of these events. In both cases, the 
i^T is typically aligned with the projection on the trans¬ 
verse plane of the second or third jet momentum. Events 
are rejected by requiring the azimuthal separation A(p 
between and (or Ei^) to be larger than 0.4. The 
resulting sample, pretag, contains 391 229 events, about 
940 of these would originate from the decay of a 300 
GeV/ W boson with SM-like couplings. 

In order to identify jets originated from the hadroniza- 
tion of a 6 quark (“6-tagged”), we use two different algo¬ 
rithms, each tuned either for making a very pure selection 
(Tight), or for making a somewhat less pure selection 
that is more efficient (Loose). The SECVTX algorithm 
m looks for a vertex displaced from the collisions point 
produced by the in-flight decay of a 6-flavored hadron; 
for this analysis we choose the tight (T) working point. 
The JETPROB algorithm m determines the probability 
that the tracks within a jet originate from the primary 
vertex; we choose for the latter algorithm the loose (L) 
working point. The efficiency for each 6-tagging algorithm 
is approximately 40-50%. We require at least one of the 
first two leading jets in Et to be tagged by the SECVTX 
algorithm. Events are further divided among twelve sta¬ 
tistically independent subsamples, depending on whether 
there are no additional 6—tagged jets (IT), or an another 
jet is tagged by jetprob but not by SECVTX (TL), or 
tagged by SECVTX (TT), the number of jets (two-jet 
or three-jet sample) and the presence or absence of a 
reconstructed electron or muon. This division increases 
sensitivity because signal-to-noise ratio and background 
composition differ across subsamples. The resulting pre¬ 
selection sample contains 25 256 events, to which a W' 
boson with SM-like couplings and 300 GeV/ (? mass would 
contribute about 480 events. 

The dominant contribution to the preselection sample 
is due to QCD multijet production. Other processes giv¬ 
ing significant contributions are top-antitop quark-pair 
production (ft), electroweak single-top-quark production, 
dibosons {WW, WZ), and production of jets in associa¬ 
tion with a boson ( V-fjets, where V stands for a IT or a 
% boson), including both heavy-flavor jets (from b or c 
quarks) and jets from light-flavor quarks or gluons that 
have been erroneously 6-tagged. 

A combination of data and simulations making use 
of Monte Carlo (MG) integration are used to derive the 
estimates for SM background contributions. The kine¬ 


matic distributions of events associated with top-quark 
pair, single top quark, T-fjets, IT -|- c, diboson iVV) 
and associated Higgs and IT or Z boson {VH) produc¬ 
tion are modeled using simulated samples. The ALPGEN 
generator [TS] is used to model T-|-jets at leading order 
(LO) in the strong-interaction coupling with up to four 
partons produced at tree level, based on generator-to- 
reconstructed-jet matching [T9l|20j. The powheg [2T] 
generator is used to model t- and s-channel single top 
quark production, while pythia [52] is used to model 
top-quark-pair, VV, and VH production. Each event 
generator uses the CTEQ5L parton distribution func¬ 
tions (PDF) [23|. Parton showering is simulated in using 
PYTHIA. Event modeling also includes simulation of 
the detector response using GEANT |23|. The simulated 
events are reconstructed and analyzed in the same way 
as the experimental data. Normalizations of the con¬ 
tributions from t- and s-channel single top quark, VV, 
VH, and tt pair production are taken from theoretical 
cross-sections |55H55] , while the normalization for IT -I- c 
production is taken from the measured cross section [29j . 
For T-Tjets production, the heavy-flavor contribution is 
normalized based on the number of 6-tagged events ob¬ 
served in an independent data control sample [30] . Con¬ 
tributions of T-fjets and VV events containing at least 
one incorrectly 6-tagged light-flavored jet are determined 
by applying to simulated events a per-event probability, 
obtained from a generic event sample containing mostly 
light-flavored jets |3ll|32]. The efhciency of the trigger- 
level selection is measured in data and applied to all 
simulated samples. 

Because QCD multijet events with large missing trans¬ 
verse energy are difficult to simulate properly, a suitable 
model is derived solely from data; we use an independent 
data sample composed of events with A(p{iBT, < 0.4 
and 50 < < 70 GeV, consisting almost entirely of 

QCD multijet contributions. First, a 6-tagging probabil¬ 
ity fi is calculated separately in each 6-tagging subsample 
i {i = IT, TL, TT) by taking the ratio between tagged 
and pretagged events as a function of several jet- and 
event-related variables [33] . Then, QCD multijet kine¬ 
matic distributions are determined separately for each 
region i by weighting the untagged data in the preselection 
sample according to the probability fi. 

The signal is modeled using pythia for IT' boson 
mass Mw' in the range 300 < Mw’ < 900 GeV/c^ in 
100 GeV je? increments, where the IT' boson is assumed 
to have purely right-handed decays. As the IT' boson 
helicity does not affect analysis observables, this model is 
valid for both a right-handed and a left-handed IT' boson 
under the assumption of no interference with SM IT boson 
production. Two scenarios are considered, depending on 
whether the leptonic decay mode IT' —>■ ^v is allowed 
or forbidden. The latter, for instance, is the case if the 
hypothetical right-handed neutrino vr is more massive 
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than the W boson. The only effect of the forbidden 
leptonic decay mode is an increased branching fraction 

b\w' tb). 

As an intermediate background-rejection step, an ar¬ 
tificial neural network, NNqcd, is employed to separate 
the dominant QCD multijet background from signal and 
other backgrounds. NNqcd is trained using event ob¬ 
servables /t[ 34]), angular observables 
A(p{^T, A^p{'^T,Ei^) and other topological informa¬ 

tions such as sphericity |3S]. As the final-state topologies 
for a W boson decaying to a top-bottom quark pair and 
s-channel single-top-quark production are similar, we em¬ 
ploy the same NNqcd function constructed to separate 
IT-l-jets events from background in the s-channel single- 
top-quark observation |36) . No information on the W 
boson mass is included in the training sample in order 
to ensure consistent performance in QCD multijet back¬ 
ground separation across the whole VD'-boson-mass range 
under study. 

The events must satisfy a minimum NNqcd require¬ 
ment to maximize sensitivity to single-top-quark s-channel 
production, which is kinematically very similar to W pro¬ 
duction at threshold. The surviving events constitute the 
signal region. To determine the appropriate normaliza¬ 
tion of QCD events in each analysis subsample, we derive 
a scale factor in the region composed by the rejected 
events. Tables|l] and |ll] show the event yields after the 
full selection. 

We use two additional neural networks, NNy-jets and 
NNj(, to classify events that satisfy the minimum require¬ 
ment on the NNqcd output variable. The first neural 
network, NNyjets, is trained to separate the W boson 
signal from D-fjets and the remaining QCD backgrounds. 
In the training, a simulated W'-boson signal is used, while 
the background sample consists of pretag data that satisfy 
the requirement on NNqcd, reweighted by the tag-rate 
probability. The second neural network, NNjf, is trained 
to separate W boson from tt production using simulated 
samples. Variables that describe the energy and momen¬ 
tum flow in the detector and angular variables are used in 
the training of the NNyjets and NN^j discriminants. The 
final discriminant, NNgig, is defined as the quadrature 
sum of the NNyjets and NN^j output variables, both mul¬ 
tiplied by an appropriate weight optimized to improve the 
expected sensitivity in each analysis subsample. Figure 
shows the expected and observed shapes of the NNgjg 
output variable for several subsamples, with the shape 
corresponding to the 300GeV/c^ W' hypothesis overlaid. 

A binned likelihood fit is performed to probe a W —>■ tb 
signal in the presence of SM backgrounds. The likelihood 
is the product of Poisson probabilities over the bins of 
the NNsig distribution. The mean number of expected 
events in each bin includes contributions from each back¬ 
ground source and from the W' —> tb process assuming a 
given value of Mw' ■ We employ a Bayesian likelihood [57] 


TABLE I: Numbers of expected and observed two-jet events 
with and without identified leptons, combined, in the IT, 
TL, and TT subsamples. The uncertainties on the expected 
numbers of events are due to the theoretical and experimental 
uncertainties on signal and background modeling. Expected 
numbers of events for a right-handed W' boson with SM-like 
couplings and a mass of 300 GeV/c^ are shown. 


Category 

IT 


TL 


TT 


s-ch. single top 

98 ± 

10 

36.4 ± 

3.8 

46.1 ± 

4.3 

t-ch. single top 

167 ± 

24 

7.3 ± 

1.1 

7.9 ± 

1.1 

tt 

457 ± 

32 

140.9 ± 

11.1 

177.4 ± 

11.7 

VV 

259 ± 

18 

28.5 ± 

2.0 

27.0 ± 

2.0 

VH 

14 ± 

1 

5.4 ± 

0.5 

7.2 ± 

0.5 

V -l-jets 

3473 ± 

901 

236.4 ± 

61.1 

156.7 ± 

38.7 

QCD 

2766 ± 

103 

220.0 ± 

16.8 

101.5 ± 

12.2 

Total background 

7235 ± 

908 

674.3 ± 

64.2 

524.5 ± 

43.0 

W (300GeV/c2) 

156 ± 

10 

59.9 ± 

4.6 

84.6 ± 

7.9 

Observed 

7128 

680 


507 



TABLE II: Same as in Table I but for three-jet events. 


Category 

IT 

TL 

TT 

s-ch. single top 

50 ± 5 

13.3 ± 1.5 

16.2 ± 1.6 

t-ch. single top 

91 ± 14 

5.8 ± 0.9 

6.9 ± 1.0 

tt 

900 ± 65 

148.2 ± 11.6 

161.6 ± 10.5 

VV 

106± 8 

9.7 ± 0.7 

7.8 ± 0.6 

VH 

6 ± 1 

1.7 ± 0.2 

2.1 ± 0.2 

V -l-jets 

1360 ± 357 

80.6 ± 21.2 

51.6 ± 13.4 

QCD 

1261 ± 64 

92.8 ± 9.4 

31.8 ± 4.6 

Total background 3774 ± 369 

352 ± 26.3 

278 ± 17.5 

W (300GeV/cQ 

80 ± 5 

23.5 ± 1.9 

28.8 ± 3.0 

Observed 

3613 

388 

274 


with a uniform, non-negative prior probability for the W 
boson production cross section times branching fraction, 
cr(pp —)■ W) X BiW' —)■ tb), and truncated Gaussian pri¬ 
ors for the uncertainties on the acceptance and shapes of 
the backgrounds. We combine the twelve signal regions 
of events characterized by different 6-tagging content, jet 
multiplicity, and presence of well-identified leptons by mul¬ 
tiplying the corresponding likelihoods and simultaneously 
taking into account the correlated uncertainties. 

Systematic uncertainties include both uncertainties on 
template normalization and uncertainties on the shape 
of the NNsig distribution. Uncertainties due to the same 
source are considered 100% correlated. These uncertain¬ 
ties apply to both signal and backgrounds, and include 
luminosity measurement (6%), 6-tagging efficiency (8 to 
16%), trigger efficiency (1 to 3%), lepton identihcation 
efficiency (2%), parton distribution functions (3%), initial- 
state and final-state simulation radiation uncertainties 
(2%) and up to 6% for the jet-energy scale [Tl]. The 
uncertainties due to finite simulation sample size, and the 
uncertainties on the normalization of the production of 
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FIG. 1: Expected and observed final discriminant distributions in the signal region. The distribution for a W' boson with 300 
GeV/ mass and SM couplings is overlaid. The signal is normalized to a cross section times branching ratio of 3pb. From left 
to right, from top to bottom: IT two-jet (a), IT three-jet (b), TL two-jet (c), TL three-jet (d), TT two-jet (e) and TT (f) 
three-jet event subsamples. 


ti (3.5%), t-channel single-top quarks (6.2%), s-channel 
single-top quarks (5%), dibosons (6%) from the theoreti¬ 
cal cross-section calculations W + c (23%) from 

the measured cross section [ITJUn]) and QCD multijet (3 
to 100%, calculated from scale factors) are not correlated. 
The production rates of events with aW or a Z boson plus 
heavy-flavor jets are associated with a 30% uncertainty. 
The shapes obtained by varying the b-tagging probability 
fi by one standard deviation from their central values 
are applied as uncertainties on the shapes of the QCD 
background. Changes in the shape of the NNgig distri¬ 
bution originating from jet-energy-scale uncertainties are 
also incorporated for processes modeled with simulations. 
An uncertainty on the 6-tagging efficiency due to differ¬ 
ent performance observed in data and simulations as a 
function of the jet Et is applied to signal distributions. 

The procedure is performed for all signal mass hypothe¬ 
ses, obtaining 95% Bayesian credibility (C.L.) upper limit 
on fT(pp —>■ W) X B{W' th) as functions of Mw, using 
the methodology described in Ref. m- The expected and 
observed upper limits are shown in Fig. The observed 
limits are compatible with the expectations calculated 
assuming that no W —>■ th signal is present in the data. 
By comparing the limits on cr(pp —>■ W) x B{W' —>■ tb) 
with the theoretical next-to-leading order calculations 


for a right-handed W boson with SM-like couplings 
we exclude W' bosons for masses less than 860 (880) 
GeV/c^ in cases where W tb decay to leptons are 
allowed (forbidden). 

For a simple s-channel-production model with effective 
coupling gw'i and assuming that couplings to light and 
heavy quarks are identical, the cross section is propor¬ 
tional to g^i. By relaxing the assumption of universal 
weak coupling, the limits on the cross section are inter¬ 
preted as upper limits on as functions of M^r'. The 
excluded region of the gw'-Mw' plane is shown in Fig. 
with gw' expressed in units of the SM weak couplings, 
5SM- For a W boson with a mass of 300 GeV/c^, the 
effective coupling is constrained at the 95% C.L. to be 
less than 10% of the W boson coupling. 

In conclusion, we perform a search for a massive reso¬ 
nance decaying to tb with the full CDF II data set, cor¬ 
responding to an integrated luminosity of 9.5 fb~^. The 
data are consistent with the background-only hypothesis, 
and upper limits are set on the production cross-section 
times branching ratio at the 95% Bayesian credibility. 
For a specific benchmark model (left-right-symmetric SM 
extension), in cases where the W —)■ t6-leptonic-decay 
mode is allowed (forbidden), we exclude W bosons with 
masses lower than 860 (880) GeV/c^. For masses smaller 
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FIG. 2: Observed and expected limits on <j{pp —> W') x 
B{W' —> tb), with ±1 ct and ±2cr credibility intervals and 
theoretical predictions for a right-handed W' boson with SM- 
like couplings in cases where the leptonic decay mode W' —>■ Iv 
is forbidden (dashed line). The CDF limits are compared with 
limits from the latest W' searches from ATLAS, CMS and 
DOUHIQ]. 


CDF Run II L = 9.5 fb ’ 



FIG. 3: Observed and expected 95% C.L. upper limits on the 
conpling strength of a right-handed W boson compared to 
the SM VF-boson coupling, Pw'/ssm, as functions of in 
cases where the leptonic decay mode W' —>■ Iv is forbidden. 
The region above each line is excluded. The CDF limits are 
compared with limits from the latest W' searches from ATLAS, 
CMS and DO |51ll0|. The vertical part in each boundary region 
of the plot represents the minimum masses for which bounds 
are quoted 

than approximately 600 GeV/c^, this search yields the 
most constraining limits to date on narrow t&-resonance 
production. 
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